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Abstract 

Shock waves generated by sharp fins, glancing across a laminar boundary layer grow- 
ing over a flat plate, are simulated numerically. Several basic issues concerning the resultant 
three-dimensional flow separation are studied. Using the same number of grid points, dif- 
ferent- grid spacings are employed to investigate the effects of grid resolution on the origin 
of the line of separation. Various shock strengths ( generated by different fin angles) are 
used to study the so-called separated and unseparated boundary layer and to establish the 
existence or absence of the secondary separation. The usual interpretations of the flow 
field from previous studies and new interpretations arising from the present simulation are 
discussed. 


I. INTRODUCTION 

In the past, fluid dynamics has been divided into two branches, theoretical and 
experimental. Continuing advances in numerical methods and in computer capabilities 
have at some point qualified computation as a separate branch of fluid dynamics. Using the 
computer as a tool, computational fluid dynamics (CFD) is able to supplement the other 
two branches and to carry out its own research, development, and further advancement in 
fluid dynamics as a field of physical science. 

A shock wave generated by one body, glancing across a boundary layer on another 
body, is a problem well suited to CFD capabilities. This is one of the most common and 
important three-dimensional (3-D) inviscid/viscous interaction problems. The problem is 
so complicated that a purely theoretical approach is almost impossible. In early studies, 
investigators have relied mainly on experimental observations of surface properties, such 
as static pressure and oil flow, to deduce the basic interaction features of the flow. With 
the advent of the supercomputer, we are now on the threshold of definitive explorations, 
finding the details of the flow-fielcl structure and the underlying physical processes. 

There are various kinds of glancing shock wave problems. Instead of dealing prob- 
lems involving complicated geometries, the present paper will focus on the problem of a 
supersonic flow over a sharp fin mounted on a flat plate, shown in Fig. 1. Recently this 
simple geometry has attracted a substantial amount of interest and has been studied ex- 
tensively, (for example, Refs. 1 - 5). However, there are still many flow-field features that 
need better explanations. 

The primary purpose of this paper is to study several basic issues regarding the 
resultant 3-D flow separation. The separation of the boundary layer on the flat plate, ( the 
sidewall), under the influence of the glancing shock wave will be specially emphasized. To 
avoid the uncertainty associated with turbulent modeling in our study, the flow is assumed 
to be laminar. The compressible Navier-Stokes equations will be solved. Cases of sharp fins 
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with various wedge angles will be investigated numerically. The results are first compared 
with the experimental and computational results of Degrez [3] for Mqq = 2.25, Re = 10 5 
and a wedge angle of 6°. Results for wedge angles of 2° and 12° will also be discussed. Using 
the same number of grid points, different grid spacings are employed to study the effects 
of grid resolution on the origin of the line of separation. Various shock strengths generated 
by different fin angles are used to study the so-called separated and unseparated boundary 
layer and the existence or absence of a secondary separation. The usual interpretations 
from previous studies and new interpretations arising from the present simulation will be 
discussed. 

It should be mentioned that this type of flow field, when turbulent, is observed to 
be unsteady (see discussion in Ref. 6). However, in the present study, the flow is assumed 
laminar and steady. The applicability of the present discussion to the turbulent cases is 
based upon the assumption that the random turbulence fluctuations are suppressed and 
the turbulent flow is steady in a “mean” sense. 

II. NUMERICAL PROCEDURE 

The governing equations of the present analysis are the time-dependent, compress- 
ible Navier-Stokes equations incorporating the concept of the thin-layer approximation 
in all three directions (Baldwin and Lomax [7]). The flow is assumed laminar and the 
wall is adiabatic. A numerical procedure developed by Hung and Kordulla [8] without 
time-splitting is used. The basic numerical scheme is MacCormack’s [9] explicit-implicit 
predictor-corrector algorithm. The solution is carried out until it converges to a steady 
state. Details of the numerical technique and boundary conditions are discussed in Ref. 8. 

Figure 2 shows a typical mesh system of 57x45x27 points for a sharp fin on a flat 
plate. The apex of the fin (at x = 0.0) is placed at a distance L — 9 cm from the flat plate 
leading edge and this distance L is used in the characteristic length in the present study. 
Here (x, y, z) and (I, J, K) are used in the conventional sense of streamwise, crossflow, 
and vertical directions. The domain of computation lies in the intervals — 1.0 < x < 3.5, 
0.0 < y < 4.5, and 0.0 < z < 1.5. Using the same number of grid points, three different 
grid spacing are employed. The coarse grid is uniform in the streamwise direction and 
geometrically stretched from the fin and plate ( in the J- and K-directions). The medium 
grid has additional geometric stretching in the streamwise direction from the apex of the 
fin. The fine grid has finer spacing, ( compared to the medium grid ), near the fin in the J- 
direction and near the apex of the fin in the streamwise direction. To avoid over-stretching 
in the outer region, several zones with different stretching factors are used in the I- and 
J-directions. Smooth transition in grid spacing is ensured from one zone to another. The 
grid spacing parameters are listed in Table 1. 

III. RESULTS AND DISCUSSIONS 

The flow to be simulated has free-stream Mach number M o© = 2.25, Reynolds num- 
ber Re = 10 5 , and free-stream temperature T ^ = 263°R. The results are at first compared 
with the experimental and computational results of Degrez for a wedge angle 9 ~ 6°. 
Figure 3 shows a comparison of surface pressures at y = 5 cm. All results are in good 
agreement, except that the computational result of Degrez shows a “dip” after the sepa- 
ration. The present results and the experimental data do not show the appearance of the 
dip in surface pressure. The fine grid result (not shown in Fig. 3) is very close to the 
medium grid result. This indicates that our grid refinement does not affect the prediction 
of surface pressure. 


2 


1 . Origin of the Line of Separation 

Reference 4 contains an extensive study of the flow-field structure of this geometry for 
turbulent boundary layer. It is well accepted that the primary separation is a consequence 
of the high pressure, recovered from the shock system, which induces flow from the hn 
surface and forces the boundary layer off the sidewall. The question is, where is the origin 
of the line of primary separation? Figure 4 shows particle traces of the result based on a 
two-equation model as described in Ref. 4 for the first mesh points above the sidewall ( 
K = 2). ( The figure is provided by Horstman [10]. ) This particle trace is constructed 
by a time integration of velocity components restricted to the plane of K = 2. Since the 
plane of K = 2 is very close to the flat plate, ( normally it would have a resolution smaller 
than the size of an oil particle), the particle integrations are treated as surface particle 
traces and are considered equivalently as a simulation of oil flow in the experiment and 
as a simulation of skin-friction lines in the theoretical approach. The “oil flow” in Fig. 
4 indicates that the line of separation originates somewhere in the plate away from the 
apex of the fin and that this feature is an open-type separation. (An open-type separation 
contains only regular points, while a closed- type separation originates from a saddle point, 
as discussed in Refs. 11 and 12.) 

Figures 5a - 5c show surface particle traces for the sequence of three grid refinements. 
For the coarse grid, the separation is an open type. As the grid spacing near the leading 
edge is refined, the starting point of the open-type separation moves and eventually the 
separation becomes a closed type. This clearly demonstrates that the grid resolution 
can affect the “calculated” topology. The coarse grid simply cannot resolve the vortex 
structure, while the fine grid can. As the wedge angle increases to 12° the vortex structure 
is large enough that the medium grid ( not shown here ) is able to reveal a closed-type 
separation. 

While the existence or not of an open-type separation is still an unanswered question, 
we believe that some of the numerically observed open-type separations (for instance, Fig. 
4) result from insufficient grid resolution. Similarly, every experiment also has resolution 
problems, such as the size of oil droplets. Some of the experimentally observed open-type 
separations may be the result of low resolution of the device and facility. 

We would like to point out that the concept of a closed 3-D separated region being 
inaccessible (see Ref. 11) is valid only in the limit of particles moving near/on body 
surface. The upstream flow particles above the surface are able to access the separated 
region behind the line of separation through the spiral nature of the separation. Indeed, 
there is no 3-D separation which is totally closed by a separation surface, as a closed 
bubble; there must be some fluid flowing in and some fluid flowing out. All 3-D separation 
surfaces are a kind of vortex sheet. 

Hereafter the fine-grid result will be used for discussion, except for cases specially 
mentioned. 


2. Secondary Separation 

The second question is the existence of a secondary separation. Experiments (for 
instance Ref. 13) very often show the appearance of another oil-accumulation line behind 
the line of primary separation. This has been interpreted as indicating the existence of a 
secondary separation, ( see Fig. 6). The plots of velocity at the first mesh point above 
the flat plate, K = 2, (Fig. 7a) also show that, in addition to the outermost primary 
separation line, (not clearly visible on the figure), there is a second “line” of clustering 
or coalescence of velocity vectors. However, based on the result of surface particle traces 
(Fig. 7b) this is not a line of separation - it is merely a demarcation between regions 
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of high and low surface skin friction. One might imagine that, in a transient stage of 
an experiment, comparatively more “oil” can be driven in by the high skin friction from 
the right and less oil carried out to the left near the region of strong variation in skin 
friction. (For convenience of discussion, here left or right refers to the orientation of one 
facing the streamwise direction.) Hence the surface may show a temporary accumulation 
of oil around this second line. Even as the wedge angle increases up to 12° (Fig. 8), our 
calculations still show no evidence of the existence of secondary separation. Indeed for the 
12° case, behind the obvious primary separation line, there is a region that the particle 
traces show strongly convergent from one side and slowly divergent away from the other 
side. ( These traces eventually converge to the primary separation line. ) Hence a high 
clustering of particle traces in that region occurs. From the plot of velocity vectors in the 
plane K = 2 (Fig. 9a), one can see that this region is associated with the drastic change 
of the surface skin friction. The flow features near the wall are strongly affected by the 
surface pressure. Their relations and connections can be seen from Fig. 9b. A strong 
surface pressure gradient induces a high skin friction and a high velocity near the sidewall. 
This results in a divergence of particle traces and leads to an appearance of the so-called 
attachment line (Fig. 9b). After the strong pressure gradient, there is a region of drastic 
change of skin friction and hence a resulting appearance of clustering of velocity vectors 
near the wall. The strong pressure gradient appears maiidy on the right of the inviscid 
shock location (Fig. 9b) and the clustering of velocity on the left (Fig. 9a). Note that the 
appearance of the clustering of particle traces (Fig. 8) does not coincide with the clustering 
of the velocity vectors near the wall. Instead, showing as a band with the clustering of 
the velocity vectors on its right, it is close to the pressure plateau region (see below). It is 
this pressure plateau region that causes the particle traces to run almost parallel to each 
other before they finally converge to the primary separation line. We would suspect that 
variations in surface-flow-visualization techniques would also result in different locations 
of temporary clustering of surface-streak lines. 

Note that Degrez’s calculation showed a noticeable ‘dip’ of pressure for the 6° case 
(see Fig. 4). When it is strong enough, this drop in pressure can significantly retard the 
primary separated flow ( passing beneath the shock system in the opposite y-direction to 
the main flow ) and lead to a secondary separation (see Fig. 12 of Ref. 3). In the present 
calculation, the pressure shows a plateau region and there is no secondary separation for 
either the 6° or 12° case. ( The 12° case has a large plateau region with a little dip of 
pressure.) As discussed above a drastic change of velocity leads to a substantial change of 
skin friction which might also result in a temporary accumulation of oil flow on the surface. 
Therefore, it is possible that an accumulation of oil flow on the surface in an experiment 
is not necessarily a line of separation contrary to the usually inference. This argument has 
also arisen in previous experiments and calculations for other geometries [14, 15]. 

A note of caution should be given here. We don’t know whether the appearance of a 
secondary oil- accumulation in an experiment is caused by a secondary separation or not. 
One possible alternative is suggested here. Furthermore, there might be other mechanisms 
in an experiment, especially for the turbulent case, that could lead to an oil accumulation 
on the surface. Further detailed and careful studies are needed to answer these questions. 

3. Absence of Separation 

Whether the boundary layer on the plate is separated or not is usually determined 
by comparing the turning angle of the limiting streamline with the glancing shock angle 
in the interaction region (Fig. 10). As shown in Figs. 5 and 8, there is clearly a line, with 
clustering of particle traces, that originates from a saddle point with turning angle greater 
than the angle of the glancing shock and the flows are separated. As the wedge angle 
decreases, one would expect that the turning angle of the skin friction line will decrease 
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and eventually become smaller than the angle of the glancing shock. The flow then will 
be classified as attached What is the change of the flowfield topology from attached to 
separated flow? This question is addressed in the following section. 

In the present paper a case with 2° wedge angle was calculated with the fine grid 
distribution. Fig. 11a shows traces of particles for which the origins are almost the same 
as those of the 6° case ( Fig. 5 ). In contrast to Figs. 5 and 8, there is no obvious 
line of convergence of particle traces and the turning of skin-friction lines is smaller than 
the glancing shock angle; this would conventionally be interpreted as an attached case. 
However, a close examination of the particle traces near the apex (Fig. lib) shows that 
actually the flow is separated. Even though it is very small, the separation also is a closed 
type, and the structure is topologically the same as that for the previous 6 and 12 c ^ s ^* 
Actually, all three cases are topologically the same as the structure of a blunt-fin flow field. 

One may imagine that, under certain conditions, this type of flow may not sepa- 
rate. However, the conventional method of interpretation using the turning angle of the 
boundary layer compared to the glancing shock angle as a criterion for separation is not 
uniformly valid, as demonstrated above. To date, the simplest and most general definition 
of a 3-D separation line, in the opinion of the author is that of Legendre [16], a line of 
separation ... has no local property. Its only characteristic is to pass through a saddle 
point” (There may be doubt as to whether or not there is a local property as pointed out 
in Refs. 17 and 18.) 

Another point should be mentioned here. Figs. 12 and 21 of Ref. 5 showed that, at 
a low wedge angle, based on the oil-flow picture the boundary layer on the sidewall was 
separated from a saddle point, as a closed type of separation, on the line of symmetiy 
near the fin apex, and then gradually became attached away from the fin. (This is not an 
open-type separation as claimed in Ref. 5.) The question arises as to where and how the 
separation ends. In the opinion of the present author the simplest explanation is that this 
is case in which the definition of separation based on turning angle and oil-accumulation 
in experimental observation fails. Because of boundary layer growth on the fin, the shock 
wave is stronger and hence the pressure rise is higher near the leading edge of the fin than 
at a position downstream. The difference in pressure rise changes the turning angle and 
degree of oil accumulation. Based on the concept of continuity, limiting streamlines would 
not join together except at a singular point. A line of separation, once it originates from a 
saddle point, will either continue going downstream to infinity or terminate at a singular 
point. In the other words, once it is separated, the flow can not gradually be reattached 
without a singular point, according to topological imperatives, 

4. Separation on Fin Surface 

As sketched in Fig. 12 of Ref. 5, even at low wedge angle there is a separation on 
the fin surface, as the high pressure flow near the fin surface tries to flow across onto the 
sidewall. The question arises as to where and how the separation starts. For a boundary 
layer to separate, an adverse pressure gradient is necessary. Plots of particle traces re- 
stricted to the plane of J = 2 and surface pressure are shown in Figs. 12a and 12b for the 
12° case. In Fig. 12a, in addition to the vortex spiral nature of the primary separation, we 
can see that the separation starts from a saddle point (a closed type). Correspondingly, 
there is a low pressure region and an adverse pressure gradient that triggers the separation 
(Fig. 12b). By overlaying Fig. 12a with Fig. 12b, one can see that the separation line is 
downstream of ( or on the righthand side of), not coincident with, the pressure minimum 
along each limiting streamline. The existence of a low pressure region can be attributed to 
the gradual decrease of total pressure in the incoming boundary layer and hence a decease 
of the pressure rise after the shock (see Fig. 12b). The small increase in pressure near the 


5 



plate results from flow stagnation. Topologically, there is even a ‘hoseshoe vortex’ (Fig. 
12a) as in the blunt-fin solution, but it is too weak to induce a significant pressure gradient. 

IV. CONCLUSION AND REMARKS 

A laminar supersonic flow over a sharp fin mounted on a flat plate has been numer- 
ically simulated. Separation of the boundary layer on the flat plate was investigated for 
various grid refinements and fin wedge angles. Several basic issues concerning 3-D steady 
separation have been discussed. 

The results of the solution have demonstrated that grid resolution can affect the 
“calculated” topology. For the coarse grid, the separation is an open type. As the grid 
spacing near the leading edge is refined, the starting point of the open-type separation 
moves and eventually the separation becomes a closed type. In the opinion of the author, 
some of the numerically observed open-type separations (for example, Fig. 4) results from 
insufficient grid resolution. Similarly, experiments have resolution problems, and some of 
the experimentally observed open-type separations may well be due to low resolution of 
the device and facility. 

Based on the computation of surface particle traces, no secondary separation is 
found in the present study. A secondary oil-accumulation line has been conjectured to be 
a demarcation between regions of high and low surface skin friction. 

In a calculation with a 2 ° wedge angle, there is no obvious line of convergent particle 
traces and the turning angles of skin-friction lines are smaller than the glancing shock angle. 
This is conventionally interpreted as an attached flow. However, a close examination of 
the particle traces near the apex has shown that actually the flow is separated, and the 
structure is the same topologically as that for the blunt-fin flow field. 
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Table 1. Three Grid Systems (57x45x27) 



coarse 

medium 

fine 

Ax min 

0.10 

0.02 

.01 

A y min 

0.0025 

0.0025 

0.0015 

A 2 min 

0.001 

0.001 

0.001 

plate leading edge 

j 1 = 4 

I = 3 

I = 3 

fin leading edge 

I = 14 

I = 22 

I = 22 
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Fig. 1 Glancing shock wave over a boundary 
layer. 



Fig. 2 Mesh system in fin and flat plate. 



Fig. 4 Surface particle traces from Horstman 
(Ref. 4) 
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Fig. 3 Comparison of surface pressure at 
y = 5 cm. 
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Fig. 5 Surface particle traces (K = 2 ) for 
a = 6°: (a) Coarse grid; (b) medium 
grid. 















Fig. 9 Velocity vectors and pressure contours 
for K = 2 (a = 6°): (a) Velocity vec- 
tors; (b) pressure contours. 



Fig. 11 Surface particle traces (K = 2) for 
a = 2°: (a) General feature; (b) near 
apex 



STREAMLINES OUTSIDE 
C T^ > BOUNDARY LAYER 



(a) ATTACHED FLOW 



(b) SEPARATED FLOW 



Fig. 12 Particle traces and pressure contours on 
the fin surface (J = 2) for a = 2°: 
(a) Particle traces; (b) pressure contours 


Fig. 10 Sketches of skin-friction line pattems- 
beneath glancing shock 
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